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Abstract 

We have investigated crystalline magnetic anisotropy in the electric field (EF) for the Fe-Pt 
surface which have a large perpendicular anisotropy, by means of the first-principles approach. 
The anisotropy is reduced linearly with respect to the inward EF, associated with the induced spin 
density around the Fe layer. Although the magnetic anisotropy energy (MAE) density reveals the 
large variation around the atoms, the intrinsic contribution to the MAE is found to mainly come 
from the Fe layer. 
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Magnetoelectric properties in the solid state are attractive for the spintronics apphca- 
tions. Through the spin-orbit interaction (SOI) the magnetic and electric properties are 
connected and consequently, the electric field (EF) allows us to manipulate the magnetic 
properties jl| and the magnetic field could control the electric polarization of materials j^. 
Such properties may be allowed in the system which looses the time-reversal and space- 
inversion symmetries, e.g., surface magnetization systems. One of the important problems 
is to control the magnetic anisotropy with the voltage or the EF. The EF induced variation 
of coercvity has been observed .n the fe.xo,„agnetie se,„,eo„ductor B a.d the large per- 
pendicular magnetization metallic surface (FePt and FePd thin layers) The strength of 
magnetic anisotropy needs in the nanoscale device so as not to loose the magnetic memory 
by the thermal fiuctuation and, meanwhile, makes difficulties in reversing the magnetization. 
The large experimental observation of EF effect on coercivity has been reported ^. The 
direct estimation of magnetic anisotropy energy (MAE) was performed for the thin film, 
Au/Fe/MgO 16|. Recent theoretical works on magnetoelectric properties in the thin iron 
films have explained the variation of MAE P, |8|, based only on a simplified Bruno's relation 
(relationship between the MAE and the atomic orbital magnetic moments) 9||. 

In order to obtain a built-up technology, the stable theoretical background for the EF 
effect is required in the basis of realistic electronic structures. The present theoretical work 
shows that the surface accumulated charge modifies the strength of magnetic anisotropy 
through the modulation of electronic state at the magnetic layer. We will discuss origins of 
the change and intrinsic contributions to the MAE. 

We have carried out first-principles electronic structure calculations [l^ which employ 
ultrasoft pseudopotentials ll| and planewave basis. Except for imposing the EF, the details 
about the method and the models are the same as in the previous study for Pt/Fe/Pt(001) 
in no EF These systems have four atomic Pt layers for the substrate and the atomic 
positions of the three bottom layers were fixed to the appropriate values of the bulk fee Pt. 
The other atoms were relaxed using the calculated atomic forces under the zero EF. All the 
atomic positions were fixed for finite EEs in MAE calculations. The atomic displacements 
induced by the EF may be a future problem for the magnetoelectric physics and the MAE 
estimation. In order to impose the EF, we have applied the scheme of effective screening 
medium (ESM) developed by Otani and Sugino 13|. In the present work, the ideal conductor 
was placed away from the Fe layer by 7.0 A. Some tiny number of electrons was added in 
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the slab for induction of the EF and at the same time induction of the counterpart charge 
at the ESM surface. The strength of EF was estimated at the front of ESM. 

The MAE was estimated from the total energy of the system with the in-plane magnetiza- 
tion ([100] direction) with respect to that of the out-of-plane (perpendicular) magnetization 
([001] direction). In the present work the MAE density, -D(r), has been introduced to 
understand the local contribution of MAE from a coarse-grain region of real space. The 
total energy E^^, where m specifies the magnetization direction of system, may formally 
be divided to the local point in real space; E^{r), where J^^iiE^{r)dr = E!^^. The atomic 
contribution in the total energy (non-local part from the pseudopotential, etc.) was rede- 
fined as a Gaussian form function centered at the atomic position. The MAE density was 
defined by -D(r) = £'[100] (r) — -Epoi] (r) and the MAE is alternatively obtained by J^^n D{r)dr. 
The atomic contributions of MAE were estimated by integrating D{r) within the atomic 
sphere with the radius of 1.3 A and the layer contributions within the layer with the plane 
boundaries determined with the midpoint of atomic coordinates along the c-direction or the 
surface normal. The z-dependent MAE density A{z) was obtained by summing D{r) up 
within the in-plane unit cell in the fixed plane normal to the surface. Further, we introduced 
another integrated MAE density; T{z) = A{z')dz' . This function helps us to capture 
the feature of A{z). The functions introduced above were calculated at a given EF (S); 
D{r,S), A{z,£), and T{z,£). 

The bulk Llo-FePt has a magnetic anisotropy along the c-direction. The total MAE 
is divided to contribution of each atom (1.96 meV and —0.01 meV in Fe and Pt atoms, 
respectively) and the interstitial region (0.67 meV/f.u.) [isl. This total MAE is also 
partitioned to the Fe and Pt layers along the c-direction, resulting in 2.83 meV and —0.21 
meV in Fe and Pt layers, respectively. These results, though the values associated with Pt 
are small and negative, have never expressed a minor contribution to magnetic anisotropy 
from Pt atoms. The intra-atomic spacial variation around Pt atoms is remarkable, reflected 
from a large spin-orbit interaction on Pt atoms and the hybridization of Pt 5d with Fe 3d 
are responsible to the large magnetic anisotropy of such systems jl6|. The spacial partition 
on MAE becomes very interesting in the surface systems and also under the EF, because the 
potential gradient against electrons is changed and the electronic structure is modulated, 
thus, which will result in the modulation of effective spin-orbit parameter [l^] . 

In Pt/Fe/Pt(001), the MAE contributes mainly from Fe layer at zero field (5.50 meV/Fe), 
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FIG. 1: Magnetic anisotropy energies (MAEs) (circles) and the number of added electrons (crosses) 
as a function of the electric field for Pt/Fe/Pt(001). The lines are obtained by the least squares fit 
to the set of data. 



amounting to 121 % of the total (5.21 meV/Fe). The main negative contribution comes from 
the capping Pt layer (—0.85 meV/Fe). Imposing the inward EF the MAE decreases, associ- 
ated with increase of the number of added electrons. The similar decrease has been reported 



in the literature 
in the Llo-FePt 



4[, where the decrease of MAE speculatively coincides with the bandfilling 
181]. Figure [1] presents the EF dependence of MAEs in Pt/Fe/Pt(001), as- 



sociated with the number of added electrons in the unit cell. The method for imposing EEs 

[isj l almost unchanged (~ O.OlmeV/Fe) the MAE at zero field, compared with the previous 

estimation with the simple periodic slab approach 12|]. The dependence of the EF on the 

number of added electrons is almost linear. This slope rate corresponds to the dielectric 

constant for vacuum (eo)', ~ l.OOeo- The decrease rate of MAE with respect to the inward 

EF is estimated to be 0.35 meV per Fe atom per V/A, corresponding to the surface MAE 

of 72 fJ/Vm. This rate is larger than the previous theoretical estimation for the metallic 

fl 

surface of iron by the factor of 3.8 (~ 19 fJ/Vm) [8|. Moreover, although the substrate of 
system is different, this slope could explain a partial contribution of the change of MAE by 
voltage in the experimental measurement Gj. 

In Figs. [2](a)-(c) the EF-induced MAE densities, A{z,Si)-A{z,0), T{z,Si)-T{z,0), and 
D{r,Si)-D{r,0), where £i=-1.97V/A, are shown along 2;-direction in Pt/Fe/Pt(001). In 
these figures, the MAE density has large strengths around Pt atoms, which reflects the large 
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FIG. 2: (color online). Electric field change of the magnetic anisotropy energy (MAE) densities 
along z-coordinate and the contour (meV/A'^) maps, (a)-(c) for Pt/Fe/Pt(001). The horizontal 
solid and dashed lines in (a)(b) indicate the atomic positions and the layer's boundaries, respec- 
tively, for better visualization. The map is shown in the (110) plane on the Fe (red ball) and Pt 
(gray ball) atoms. 



SOI at them. However, such contributions to the MAE seem to cancel out around Pt atoms 
and, as shown in Fig. Mjo), the intrinsic contribution to the total induced variation on MAE 
is raised around the Fe layer. This indicates that the local electronic structure centered at 
the Fe layer is intrinsically important for the change of MAE. 

To see the relationship between the variation of MAE and the EF in more details, we 
have calculated the induced spin density, the energy-dependent density of states for five 3d 
angular-dependent orbitals on Fe atom, and the orbital-specified band dispersions [l2 1 . Fig- 
ures [3] (a) and (b) show the induced majority- and minority-spin densities along ^-direction 
and the induced spin density map, respectively, at S = — 1.97V/A in Pt/Fe/Pt(001). At 
Fe the minority-spin density is induced and this induction is found to result in a partial 
cancellation of the increase of (i3^2_^2 and dx2_y2 components against the decrease of d^z 
and dyz. The latter is associated with the charge accumulation at outside of the capping 
layer through the orbital hybridization between Fe 3d and Pt 5d, which implies a subtle 
bonding reduction. As a result, the d^z and dy^ components grow just above the Fermi level. 
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FIG. 3: (color online). The induced majority- (solid curves) and minority- (dashed curves) spin 
charge densities along z-coordinate and the contour maps of induced spin densities (xlO"'^e/A^) 
(a)(b) for Pt/Fe/Pt(001). The magnetization direction used for the figure drawing is perpendicular 
to the surface plane. See the caption at Fig. [2] for other explanations. 



These EF modulations in electronic structure around the Fermi level can be related with 
the decrease of MAE by using the second order perturbative consideration loj l; 

MAE oc ^^ KM4|k^)|^-|(ko|4|W|^ ^ 

k o,u ^kw ^ko 

where ko and ku specify the occupied and unoccupied states with the wave vector k and 
£o(a = x,z) the angular momentum operators. The EF modulation mentioned above in- 
creases the couplings between occupied and unoccupied states (couplings of {Sz"^ — r'^\\yz) 



and {yz\\x'^ — y'^)) through the operator and, thus, reduces the MAE |19|. 

In Table [U the spin and orbital moments are reported. These quantities change linearly 
with the EF in Pt/Fe/Pt(001). One would explain the MAE from the orbital magnetic 
moments, while it is not probable that the simplified Bruno's relation applies to the FePt 

n n 

alloy (see Eq.(9) in [IJ]) without the spin flip contribution [9|]. This is because the Ft 
does not have any large exchange splitting. As implied by the difference of atomic orbital 
magnetic moments on Fe and Pt atoms in Table [U the MAE may be supposed to mainly 
come from the Pt atoms. This picture is in contradiction with the feature obtained from the 
MAE density (Fig. [2]). In this context, to give a reasonable explanation in the relation with 
the orbital magnetic moment and the MAE, the application of a general Bruno's relation 
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is required, accompanied with the spin flip contribution and the interstitial (inter-atomic) 
contribution to the MAE. 

Imposing the EF, the number of electrons changes in a few layers of the surface due to 
a screening effect of metal, as shown Fig. 3(a). This feature should be realized also in the 



experiment 



a. 



If the change of coercivity is assumed to be attributed to the one or two 
magnetic layers of the metal surface in the experiment, the MAE of these layers could reduce 
by 10 ~ 30 % for the EF of — 0.03V/A. This reference to the experiment provides a good 
agreement with our decrease of MAE (13 % in Pt/Fe/Pt(001)) by the EF of -0.030V/A 
which is scaled by the dielectric constant of the experimental substrate on metal surface 
(— 1.97V/A on vacuum). 

In summary, we have studied the EF dependence of MAE for the large perpendicular 
magnetic anisotropy surface system. The MAE linearly decreases with the inward EF in 
Pt/Fe/Pt(001). By analysis of the induced MAE density, the large variation around Pt 
atoms was revealed and the intrinsic contribution to the MAE was verified to mainly come 
from the Fe layer. The analysis of MAE density provides a promising tool for pursuing the 

TABLE I: Magnetoelectric properties in Pt/Fe/Pt(001); magnetic anisotropy energies(MAEs) at 
zero field and their gradients with respect to the electric field(i-), the spin magnetic moment (SMM) 
from the three surface layers in the [001] magnetization, and the difference of atomic orbitals 
magnetic moments(DOMM) between the [001] and [100] magnetizations, Morb[001] - Morb[100]. 



Pt/Fe/Pt(001) 


MAE at £ = 


5.21 (11)" 


MAE slope rate 


0.35 (72)^ 


SMM at ^ = 


3.72 " 


SMM slope rate 


0.0065(0.99) 


Fe DOMM at £: = 


0.0080 


Fe DOMM slope rate 


-0.0037 (-0.56) 


Pt(c) DOMM at £: = 


-0.0417 ^ 


Pt(c) DOMM slope rate 


-0.0024 (-0.36) 



" in meV/Fe (J/m^), in meV/Fe per V/A (fJ/Vm), 
" in ^B, in per Y/k{10~^^ GmVV) 
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spacial contribution of nanoscale structures, regardless of the knowledge of orbital magnetic 
moments. The present study indicates that the relative modification in the electron filling 
of each 3d orbital by the EF, resulting in the accumulated charge at the magnetic layer, 
causes the variation of MAE. 
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